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We analyze the physics of spin-3 Bose-Einstein condensates, and in particular the new physics 
expected in on-going experiments with condensates of Chromium atoms. We first discuss the ground- 
state properties, which, depending on still unknown Chromium parameters, and for low magnetic 
fields can present various types of phases. We also discuss the spinor-dynamics in Chromium spinor 
condensates, which present significant qualitative differences when compared to other spinor con- 
densates. In particular, dipole-induced spin relaxation may lead for low magnetic fields to transfer 
of spin into angular momentum similar to the well-known Einstein-de Haas effect. Additionally, a 
rapid large transference of population between distant magnetic states becomes also possible. 

PACS numbers: 



Within the very active field of ultra cold atomic 
gases, multicomponent (spinor) Bose-Einstein conden- 
sates (BECs) have recently attracted a rapidly growing 
attention. Numerous works have addressed the rich va- 
riety of phenomena revealed by spinor BEC, in partic- 
ular in what concerns ground-state and spin dynamics 
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The first experiments on spinor BEC were performed 
at JILA using mixtures of 87 Rb BEC in two magneti- 
cally confined internal states (spin-1/2 BEC) .10]. Op- 
tical trapping of spinor condensates was first realized in 
spin-1 Sodium BEC at MIT constituting a major 
breakthrough since, whereas magnetic trapping confines 
the BEC to weak-seeking magnetic states, an optical trap 
enables confinement of all magnetic substates. In addi- 
tion, the atoms in a magnetic substate can be converted 
into atoms in other substates through interatomic inter- 
actions. Hence, these experiments paved the way to- 
wards the above mentioned fascinating phenomenology 
originating from the spin degree of freedom. Various ex- 
periments have been realized since then in spin-1 BEC, 
in particular in 87 Rb in the F = 1 manifold. It has 
been predicted that spin-1 BEC can present just two dif- 
ferent ground states phases, either ferromagnetic or po- 
lar 0. In the case of F = 1 87 Rb it has been shown 
that the ground-state presents a ferromagnetic behav- 
ior 0, Q, 0| . These analyses have been recently ex- 
tended to spin-2 BEC, which presents an even richer va- 
riety of possible ground-states, including in addition the 
so-called cyclic phases 0, 0- Recent experiments have 
shown a behavior compatible with a polar ground-state, 
although in the very vicinity of the cyclic phase 0, ^| . 
Recently, spin dynamics has attracted a major interest, 
revealing the fascinating physics of the coherent oscilla- 
tions between the different components of the manifold 
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Very recently, a Chromium BEC (Cr-BEC) has been 



achieved at Stuttgart University 20]. Cr-BEC presents 
fascinating new features when compared to other exper- 
iments in BEC. On one hand, since the ground state of 
52 Cr is 7 S3, Cr-BEC constitutes the first accessible exam- 
ple of a spin-3 BEC. We show below that this fact may 
have very important consequences for both the ground- 
state and the dynamics of spinor Cr-BEC. On the other 
hand, when aligned into the state with magnetic quan- 
tum number m — ±3, 52 Cr presents a magnetic moment 
/i = 6/is, where /x_b is the Bohr magneton. This dipo- 
lar moment should be compared to alkali atoms, which 
have a maximum magnetic moment of \\ib , and hence 36 
times smaller dipole-dipole interactions. Ultra cold dipo- 
lar gases have attracted a rapidly growing attention, in 
particular in what concerns its stability and excitations 
|21j . The interplay of the dipole-dipole interaction and 
spinor-BEC physics has been also considered 0, • R- e ~ 
cently, the dipolar effects have been observed for the first 
time in the expansion of a Cr-BEC |22^ . 

This Letter analyzes spin-3 BEC, and in particular 
the new physics expected in on-going experiments in Cr- 
BEC. After deriving the equations that describe this sys- 
tem, we focus on the ground-state, using single-mode 
approximation (SMA), showing that various phases are 
possible, depending on the applied magnetic field, and 
the (still unknown) value of the s-wave scattering length 
for the channel of total spin zero. This phase diagram 
presents certain differences with respect to the diagram 
first worked out recently by Diener and Ho [2^]. In the 
second part of this Letter, we discuss the spinor dynam- 
ics, departing from the SMA. The double nature of Cr- 
BEC as a spin-3 BEC and a dipolar BEC is shown to 
lead to significant qualitative differences when compared 
to other spinor BECs. The larger spin can allow for fast 
population transfer from m = to m = ±3 without a 
sequential dynamics as in F — 2 87 Rb In addition, 
dipolar relaxation violates spin conservation, leading to 
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rotation of the different components, resembling the well- 
known Einstein-de Haas (EH) effect [3 . 

In the following we consider an optically trapped Cr- 
BEC with N particles. The Hamiltonian regulating Cr- 
BEC is of the form H = H + V sr + V dd . The single- 
particle part of the Hamiltonian, H , includes the trap- 
ping energy and the linear Zeeman effect (quadratic Zee- 
man effect is absent in Cr-BEC), being of the form 
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H = J dvJ2^l(r) 
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i>m (r), 
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where ip^ (ipm) is the creation (annihilation) operator in 
the m state, M is the atomic mass, and p = g/isB, with 
g = 2 for 52 Cr, and B is the applied magnetic field. 

The short-range (van der Waals) interactions are given 
by V sr - For any interacting pair with spins Si j2 , V sr 
conserves the total spin, S, and is thus described in terms 
of the projector operators on different total spins Vs, 
where S = 0, 2, 4, and 6 (only even S is allowed) Q : 



If 6 
V sr = r / drJ2gsPs(r) 
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where gs = 'inh 2 a s / 'M , and as is the s-wave scattering 
length for a total spin S. Since Si • S 2 = (S 2 - Sf - S%)/2, 



then Es^s(r) 



! (r) :, Es X sV s (r) =: F 2 (r) :, and 



J2g Ag'Ps(r) =: 2 (r) :, where :: denotes normal order, 
A S = [S(S + 1) - 24J/2, n(r) = ^(r)^(r), F 2 = 

£i=W^?« with M T ) = E ro „^( r )^mn^n( r )> and 

O 2 = V , (>: r with Oij = J2 m ,n^L(r)(S l S% m Mr)- 
gx,y,z are ^ e g pj n ma trices. Hence, 



dr[c : n 2 (r)+c 1 F 2 {r) + c 2 Vo(r)+c 3 6 2 {r) :], 

(3) 

where V Q {r) = \ ^(-rt^W-., and c = 
(-llfla + 8Ig 4 + 7 56 )/77, ci = (<? 6 - g 2 )/18, c 2 = 
5o+(-55 52 +27 ff 4-5 5 6)/33, c 3 = ff2 /126- ff4 /77+g 6 /198. 
For the case of 52 Cr |25j, &6 = 112as, where as is 
the Bohr radius, and cq = 0.65<?6, ci = 0.059(76, c 2 = 
,g0 + 0.374^6, and C3 = — 0.002<?6- The value of ao is 
unknown, and hence, in the following, we consider differ- 
ent scenarios depending on the value of go/ge- Note that 
Eq. J3J) is similar to that obtained for spin-2 BEC 0, , 
the main new feature being the C3 term, which introduce 
qualitatively new physics as discussed below. 
The dipolc-dipole interaction Vdd is of the form 



i;,v, - | / dr j dr' T — ^^(r)^(r') 
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- 3(S„ m • e)(S m , n , • e)] ^(r)^(r'), (4) 

where Cd — iiaii%g s /kK, with /j,q the magnetic perme- 
ability of vacuum, and e = (r — r')/|r — r'|. For 52 Cr 
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FIG. 1: Phase diagram function of go/ge and p/ge, where 
p = 2gsn B B/NI3. 



Cd = 0.004<76- Vdd may be re-written as: 

dvdt 
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r — r 

^ 2 ,_(r,r')r 2 i+^, + (r,r')r 2 _i 
^_,_(r,r')r 22 +^ + . + (r,r')F 2 - 2 } :, (5) 



where F zz {r,r') = ^2/3{3F z (r)F z (r') - F(r) • 
F(r')], ^,±(r,r') = ±[F ± (v)F z (v>) + F z (v)F ± (v% 

^±,±(r. r') = F±(r)F±(r'), F± = F x ± %F y , and Y 2m (r - 
r') are the spherical harmonics. Note that contrary to 
the short-range interactions, Vdd does not conserve the 
total spin, and may induce a transference of angular mo- 
mentum into the center of mass (CM) degrees of freedom. 

We first discuss the ground-state of the spin-3 BEC 
for different values of go, and the magnetic field, p. 
We consider mean-field (MF) approximation ^ m (r) ~ 
^/Nip m (r). In order to simplify the analysis of the possi- 
ble ground-state solutions we perform SMA: VVn( r ) = 
<f>(r)^ m , with J<ir|<I>(r)| 2 = 1. Apart from spin- 
independent terms the energy per particle is of the form: 
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- 2r + i/ 2 /_ + 2r_i/ 2 / + + r +2 / 2 +r_ 2 / 2 ] , (6) 

where = J dr|$(r)| 4 , S _ = |E ra (-l) m ^ m i m , /, = 

Em,n Tpm Si mni>n, O tj = J2 m ,n ftf^ Wn and T m = 

/ dr / dr' I $ (r') 1 2 1 $ (r ) 1 2 r 2m (r - r' ) / 1 r - r' 1 3 . Note that 
r±i = for any symmetric density |<&(r)| 2 . 

Let us consider a magnetic field in the z-direction. 
Hence, the ground-state magnetization must be aligned 
with the external field, and /+ = /-= 0. Then, 
we obtain that O? = 77 - 12 72 + 37 2 /2 - 
/ 2 /2 + M 2 /2 + 2|q| 2 , where lz = E m m 2 |V> m | 2 , r) = 
EJU-3 VI12 - (m + 2)(m + 1)][12 - m(m + l)]^ +2 Vm 
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moving spin- independent terms, the energy becomes 
E = N0e/2, with 
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where p = 2p/N0, and Ci = C\ — c 3 /2 + ^/167r/5r Cd//3. 
Since c<j <C c l5 dipolar effects are not relevant for the 
equilibrium discussion. We will hence set To = 0. 

We have minimized Eq. JJJ) with respect to ip m , un- 
der the constraints J2 m IV'ml 2 = 1 an d /+ = |26| . 
Figs. H shows the corresponding phase diagram, which 
although in basic agreement with that worked out re- 
cently in Ref. j 2 ■' -J | presents certain differences in its final 
form. For the phases discussed below a = 27]. For 
sufficiently negative 50 the system is in a polar phase 
P = (c6>, 0, 0, 0, 0, 0, s0), with c = cos, and s = sin. This 
phase is characterized by f z ~ —3p/p c (p c ~ 65i ~ 0.36), 
4|s_| 2 = 1 - (p/p c ) 2 , lz = 9, and \ v \ = 0. The po- 
lar phase extends up to 50/56 = 0.01 for p = 0. For 
sufficiently large 50 an d p, CY-ap = (cd, 0, 0, 0, 0, s9, 0) 
occurs. This phase is cyclic (s_ = 0) and it is charac- 
terized by f z = -3p/p c , 7 Z ~ 6 + 3p/p c , and \r/\ = 0. 
Both P and Cy_ 3 . 2 continuously transform at p = p c 
into a ferromagnetic phase F = (1,0,0,0,0,0,0). For 
0.04 < p < 25i + 9c 3 CF_ 3i2 becomes degenerated with 
the cyclic phase CY_i )2 = (0, 0, c9, 0, 0, 0, s0, 0). The 
latter state differs from since 7 Z = 2 — 3p/p c . 

For sufficiently large go and p < 0.04, a cyclic phase 
(CY_3_i,i )3 ) of the form (i/>- 3 , 0, V'-i, 0, 0, V-3) oc- 
curs. Contrary to the other cyclic phases this phase is 
characterized by \r]\ > |28(. Finally for a region around 
go = an additional phase with \r)\ > 0, |s_| > is 
found. In this phase two possible ground states are de- 
generated [13, namely 5_ 2 ,o, 2 = (0,-0-2, 0,^0,0,^+2,0), 
and -1^1.3, which has a similar form as C7Y1-3 —1^3. 

For a Cr-BEC in a spherical trap of frequency lo in 
the Thomas-Fermi regime, p = p c for a magnetic field (in 



mG) B ~ 1.25 (TV/10 5 ) /0 (cj/10 3 x 2vr) /j . Hence, most 
probably, magnetic shielding seems necessary for the ob- 
servation of non-ferromagnetic ground-state phases. We 
would like to stress as well, that there are still uncertain- 
ties in the exact values of 02.4,6- In this sense, we have 
checked for a§ — 98as, = 6Aas, and a 2 = — 27as, that 
apart from a shift of 0.2 towards larger values of 50/56, 
the phase diagram remains qualitatively the same. 

In the last part of this Letter, we consider the spinor 
dynamics within the MF approximation, but abandoning 
the SMA. The equations for the dynamics are obtained 
by deriving the MF Hamiltonian with respect to i\)* m (r) : 
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t- r± 2 , T _ 

and S, 



where .Ao = y/6n/5[y/8/3To, z -r .1 1 . 1- 1 

At = v / W5[- r o,±A/6 t r ±1 , 

J dr7 i (r / )y 2m (r - r')/|r - r'| 3 , 

yT2 — m(m =F 1)< Note that all if} m , n, ft, Oij, and A% 
have now a spatial dependence. We first consider p = 0, 
and discuss onp 7^ below. There are two main features 
in the spinor dynamics in Cr-BEC which are absent (or 
negligible) in other spinor BECs. On one side the c 3 
term allows for jumps in the spin manifold larger than 
one, and hence for a rapid dynamics from e.g. m = to 
m = ±3. On the other side, the dipolar terms induce a 
EH-like transfer of spin into CM angular momentum. 
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FIG. 2: |i/; m= _ 2 (r)| 2 at ut = 40 (a) and 120 (b) for p = 0, 
go = 0, Lu z = 1kHz, N = 10 4 atoms, and ip(t = 0) = ipm=-3- 
The x and y axes are in y^h/muj units. 

We consider for simplicity a quasi-2D BEC, i.e. a 
strong confinement in the z-direction by a harmonic po- 
tential of frequency u> z . Hence if> m (r) = (f>o(z)ip m (p), 
where 4>o{ z ) — exp[—z 2 /2l 2 ]/n 1 / 4 y/T z ~, with l z = 
yJhJmJuT z . We have then solved the 2D equations us- 
ing Crank-Nicholson method, considering a harmonic 
confinement of frequency u> in the xy-plane. In 2D, 
r±i z = 0, but these terms vanish also in 3D due to 
symmetry, and hence the 2D physics is representative of 
the 3D one. The vanishing of T±i iZ is rather important, 
since if ip(t — 0) = ip m =±3, ^±i.z is responsible for a 
fast dipolar relaxation (for u>t ~ 1). Hence, a BEC with 
ip(t = 0) = ip±3 does not present any significant fast spin 
dynamics. However significant spin relaxation appears in 
the long time scale, due to the T 2j _(r) terms. This is the 
case of Fig. |2 where we ip(t = 0) = ipm=-3 |3(ij ]. 

One of the most striking effects related with spin re- 
laxation is the transference of spin into CM an gula r mo- 
mentum, which resembles the famous EH effect [24j . The 
analysis of this effect motivated us to avoid the SMA in 
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FIG. 3: Population of V'0,1,2,3 versus uit for go — ge, p = 0, 
and V'(O) = tjjQ. Note the rapid growth of m = ±3 (arrow). 

the study of the spinor dynamics. In Figs. El we show 
snapshots of the spatial distribution of the m = — 2 com- 
ponent. Observe that the wavefunction clearly loose its 
polar symmetry, since spin is converted into orbital an- 
gular momentum. The spatial patterns become progres- 
sively more complicated in time. 

The other special feature of Cr-BEC, namely the ap- 
pearance of the C3 term can have significant qualitative 
effects in the dynamics both for short and for long time 
scales. The evolution at long time scale may present in- 
teresting features, as large revivals, and it will be consid- 
ered in future work. Here, we would like to focus on the 
short-time scales, where the C3 term may produce fast 
transference from ip m =o to the extremes ip m —±.^. The 
latter is illustrated in Fig. EI where we consider go = g$. 
The population is initially all in the m = [3(| ■ Contrary 
to the case of F = 2 in 87 Rb |l3|, there is at short time 
scales a jump to the extremes (the population of ±m is 
the same due to symmetry since p = 0). This large jump 
is absent if C3 = 0, and depends on the value of go. In 
particular, if go — one obtains at short time scales a 
sequential population as for F = 2 87 Rb fl3j |. 

We finally comment on the dynamics if p 7^ 0. In the 
case of F — 1 or F = 2 87 Rb, the dynamics is indepen- 
dent of p since the linear Zeeman effect may be gauged 
out by transforming ip m — > ip m e lpmt / n , due to the con- 
servation of the total spin. In Cr-BEC the situation is 
very different, since the r±i jZ and F±2,=f do n °t conserve 
the total spin, and hence oscillate with the Larmor fre- 
quency ujl = p/h and 2u>l, respectively. If u>l 3> w one 
may perform rotating-wave approximation and eliminate 
these terms. Hence the coherent EH-like effect disappears 
for sufficiently large applied magnetic fields. 

In conclusion, spin-3 Cr-BEC is predicted to show dif- 
ferent types of spin phases depending on ao and the mag- 
netic field. The spinor dynamics also presents novel fea- 
tures, as a fast transference between tpo — ► ^±3; an d the 
Einstein-de Haas-like transformation of spin into rotation 
of the different components due to the dipole interaction. 
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